Source for Intra-Pulse Multi-Energy X-Ray Cargo Inspection

ABSTRACT

Methods for generating a multiple-energy X-ray pulse. A beam of electrons is generated with an electron gun and modulated prior to injection into an accelerating structure to achieve at least a first and specified beam current amplitude over the course of respective beam current temporal profiles. A radio frequency field is applied to the accelerating structure with a specified RF field amplitude and a specified RF temporal profile. The first and second specified beam current amplitudes are injected serially, each after a specified delay, in such a manner as to achieve at least two distinct endpoint energies of electrons accelerated within the accelerating structure during a course of a single RF-pulse. The beam of electrons is accelerated by the radio frequency field within the accelerating structure to produce accelerated electrons which impinge upon a target for generating Bremsstrahlung X-rays.

The present application is a continuation application of U.S. patent application Ser. No. 15/832,369, entitled “Source for Intra-Pulse Multi-Energy X-Ray Cargo Inspection” and filed on Dec. 5, 2017, which is a continuation application of U.S. patent application Ser. No. 15/307,463, of the same title, filed on Oct. 28, 2016, and issued as U.S. Pat. No. 9,867,271 on Jan. 9, 2018, which is a 371 National Stage Entry of PCT Application No. PCT/US15/30716, of the same title and filed on May 14, 2015, which in turn, claims the priority of U.S. Provisional Application Ser. No. 61/994,484, of the same title and filed on May 16, 2014. All of the above applications are incorporated herein by reference in their entirety.

TECHNICAL FIELD

The present invention relates to systems and methods for inspection of cargo using penetrating radiation, and, more particularly to systems and methods for material discrimination on the basis of varying both the energy and the flux of incident radiation during the course of a single pulse.

BACKGROUND OF THE INVENTION

Discrimination of materials has become a standard requirement for security inspection systems. An inspection system for cargo and containers screening typically employs an electron accelerator capable of interlaced dual energy operation, and differential transmission of X-rays characterized by distinct energy spectra can be used to distinguish among intervening materials of different atomic numbers. The term “interlaced energies,” as employed herein, denotes the use of a stream of X-ray pulses wherein successive pulses are characterized by distinct energy spectra. The use of interlaced energy inspection for material discrimination is well known, with processing techniques surveyed, for example, by Ogorodnikov et al., Processing of interlaced images in 4-10 MeV dual energy customs system for material recognition, Phys. Rev. Special Topics—Accelerators and Beams, vol. 5, 104701 (2002), and in references cited therein, all of which are incorporated herein by reference. A Bremsstrahlung spectrum is characterized by its endpoint energy, defined by the energy of electrons impinging upon an X-ray target in order to generate the X-rays. Attenuation by an inspected object of transmitted radiation for two (or, on some occasions, more) distinct energies provides the input data for identifying the type of material being inspected.

Various techniques are known for generating X-rays of interlaced energies based on electron accelerators, such as described, for example, in U.S. Pat. No. 7,646,851, entitled “Device and Method for generating X-Rays Having Different Energy Levels and Material Discrimination System,” and U.S. Pat. No. 8,604,723, entitled “Interlaced multi-energy radiation sources,” both of which are incorporated herein by reference. The technologies of interlaced energy irradiation merit no further discussed here, since they are irrelevant to the present invention, other than to highlight their deficiencies for purposes of cargo inspection.

Several limitations accompany material discrimination functionality when a source of interlaced X-ray energy is employed. Using two pulses separated in time for producing one inspection data point effectively reduces inspection speed. Moreover, while the basic assumption of dual-energy techniques is that the same region of the cargo is being probed by both energies, it must be borne in mind that the cargo and the probe are typically in relative motion. Interlaced energy approaches are thus only tenable for slow moving cargo.

X-ray security inspection systems for the inspection of cargo and shipping containers typically use transmission radiographic techniques. FIG. 1 depicts a cargo inspection system employing such a technique. A fan-shaped beam 12 of penetrating radiation, emitted by a source 14, is detected by elements 18 of a detector array 16 distal to a target object, here truck 10, is used to produce images of the target object. Detector elements 18 produce corresponding detector signals processed by processor 19 to provide information relative to the material composition of the cargo and images of its spatial distribution. The thickness of material to be penetrated by the X-rays may exceed 300 mm of steel equivalent in some cases. To insure the required penetration, inspection systems typically use X-rays with a maximum energy of several MeV, currently up to about 9 MeV. X-rays in excess of 1 MeV are frequently referred to as hard X-rays or high-energy X-rays. While the invention described herein pertains to any penetrating radiation, it may be described, purely as a matter of heuristic convenience, in terms of high-energy X-rays.

Information (such as mass absorption coefficient, effective atomic number Z_(eff), electron density, or the spatial distribution of any of the foregoing, etc.) with respect to the material composition of the contents of objects may be obtained on the basis of the interaction of X-rays with the material, and, more particularly, by illuminating the material with X-ray beams having energy spectra with more than one distinct energy endpoint (peak energy), or by employing energy discriminating detectors. Dual energy methods of material discrimination are widely used in X-ray inspection systems for security control of hand luggage in customs and other security checkpoints.

Dual (and, more generally, multiple-) energy methods have been extended to high-energy inspection systems for cargo containers, where material discrimination is less effective due to the weaker Z-dependence of the dominant interaction.

In the practice of dual-energy inspection, X-ray transmission data of an inspected object are obtained for both energies, and processed by computer, whereupon a resulting image is displayed on a monitor, typically in a special color palette that facilitates visual identification of contraband or hazardous materials. More particularly, special computer software may identify various materials and artificial colors may be assigned to various values of Z_(eff).

A typical energy range for the inspection of smaller objects is below 0.5 MeV, taking advantage of the strong Z-dependence of the X-ray attenuation coefficient due to the prevalence of the photoelectric interaction (characterized by a cross-section, ˜Z⁴-Z⁵) at lower energies. In the range of 1-10 MeV, however, X-ray interaction is dominated by the Compton effect with its weak dependence of attenuation coefficient (mass absorption) on the atomic number: μ_(c)˜Z/A (which is approximately constant and equal to 0.5), where Z denotes atomic number, and A denotes atomic mass, which is to say that the mass absorption coefficient is largely Z-insensitive in the energy regime dominated by Compton scatter.

A preferred method for material discrimination entails variation of the pulse energy during the course of each single pulse, as described in detail in U.S. Pat. No. 8,457,274 (“Arodzero '274”, issued Jun. 4, 2013), which is incorporated herein by reference.

Leó Szilárd conceived of the linear accelerator (linac) in 1928, while a professor at the University of Berlin. A linac was also constructed independently by Rolf Widerøe, then an engineering graduate-student under Walter Rogowski at Aachen, at about the same time. Electrons accelerated by a linear accelerator were first used to generate X-rays at Stanford in the mid-1950's.

Some prior art methods for varying the emitted energy during the course of a pulse have required that the x-ray flux track the end-point energy. The Arodzero '274 Patent, for example, states that “Concurrently with the sweeping of the endpoint energy, the X-ray flux may increase from a minimum to a maximum.” (Arodzero '274, col. 6, lines 47-48.)

US Published Patent Application 2014/0270086 (to Krasnykh), incorporated herein by reference, describes an intra-pulse multi-energy method that uses a traveling wave accelerator structure. It suggests the use of feedback to the electron gun grid voltage to compensate for X-ray flux variation during the course of a pulse. Krasnykh et al., Concept of RF Linac for Intra-Pulse Multi-Energy Scan, SLAC Pub-15943, (Apr. 18, 2014) provides further description, and is also incorporated herein by reference. The prior art mode of operation, however, could not accommodate separate tailoring of the flux and end-point energy of an X-ray pulse, even though such operation would be highly advantageous in a cargo inspection context, for example.

One of the limiting factors of inspection speed is RF-power available for accelerating. The maximum pulse repetition frequency (PRF) that a linac-based X-ray source can provide is limited by the RF source. The RF source (typically, a magnetron or a klystron) has limitations on maximum average P_(av,max) and pulsed P_(p,max) power. These two parameters define the maximum duty factor dmax, which also can be expressed in terms of PRF (f) and pulse duration t_(p):

$\begin{matrix} {d_{m\; {ax}} = {\frac{P_{{av},{{ma}\; x}}}{P_{{p,{{ma}\; x}}\;}} = {f \cdot {t_{p}.}}}} & (1) \end{matrix}$

For example, where a single energy (SE) (non-interlaced) accelerator, characterized by P_(p,max) is chosen to produce the high energy (HE) pulse, with t_(p)≈3.3 μs and d_(max)≈0.001, the maximum PRF would be limited to f_(H)≈300 Hz (pps).

For a dual-energy interlaced linac, the maximum available frequency can be estimated from the equation

$\begin{matrix} {{{f_{DE} \approx \frac{P_{{av},{{ma}\; x}}}{P_{H} \cdot t_{p} \cdot \left( {1 + {P_{L}/P_{H}}} \right)}} = \frac{f_{H}}{\left( {1 + {P_{L}/P_{H}}} \right)}},} & (2) \end{matrix}$

where P_(H) and P_(L) represent the RF power necessary to produce high (HE) and low (LE) energy pulses, respectively. If the assumption is made that t_(p) remains the same for both energies, and that P_(H)=P_(p,max), then, for P_(L)=P_(H) (RF-power remaining constant for both pulses, achieved, for example, by RF-switches/regulators, manipulation of beam loading, and phase-shifting of the accelerating field), f_(DE)=1/2·f_(H). That is to say, that a dual energy repetition rate of, at best, half that of the single-energy rate, may be achieved. On the other hand, if the low-energy pulses produce only half the power of the high-energy pulses, P_(L)=1/2PH, (as might be implemented using RF-generator power supply modulation, for example), then f_(DE)=2/3 ·f_(H), which is to say that ⅔ of the single-energy pulse rate may be achieved on an interlaced energy basis.

In prior practice, both the RF-power and the injected beam were turned on at the same time (t_(b)=0). The result of such prior art practice is shown by the dotted curve 30 plotted in FIG. 3, which represents the beam energy W vs. time t for a 6-MeV accelerating structure designed for security applications. The filling time, which is the time it takes for the electric field in the accelerator structure to decay to e⁻¹ of its initial value, is t_(f,95%)≈1 μs .

A well-known technical solution for reducing the filling time was described by Roger Miller, Comparison of Standing-Wave and Travelling-Wave Structures, SLAC Linear Accelerator Conference, SLAC-PUB—3935, pp. 216-21 (1986) (hereinafter, “Miller (1986)”, which is incorporated herein by reference. The Miller solution allows for creating the beam pulse with constant energy over the pulse duration. The accelerating beam turns on with delay t_(b) that is defined as:

$\begin{matrix} {t_{b} = {\tau \cdot {{\ln\left( \frac{\sqrt{4\beta \; {rLP}}}{IrL} \right)}.}}} & (3) \end{matrix}$

βis the coupling coefficient between an RF power feed waveguide 222 (shown in FIG. 2) and an accelerating structure 22 (shown in FIG. 2, also referred to herein as an “accelerating cavity structure”), r is the shunt impedance of the accelerating structure 22, L is the length of the accelerating structure 22, and P is the power dissipated in the accelerating structure 22, and τ is the decay time constant of the accelerating structure 22. (Thus, both numerator and denominator of the logarithmic argument have units of voltage.) Rigorously, is defined as the ratio of power lost outside the accelerating cavity structure 22 (i.e., in the feed waveguide 222) to the power dissipated inside the accelerating cavity structure 22. If β=/β₀ has been adjusted so that there is no RF-power reflection from the accelerating structure 22 when the beam 220 is on, the above equation can be cast as:

$\begin{matrix} {t_{b} = {\tau \cdot {\ln \left( \frac{2 \cdot \beta_{0}}{\beta_{0} - 1} \right)}}} & (4) \end{matrix}$

where β₀ is the optimum coupling coefficient,

$\begin{matrix} {{\beta_{0} = \left( {{\frac{I}{2}\sqrt{\frac{rL}{P}}} + \sqrt{1 + {\frac{I^{2}}{4}\frac{rL}{P}}}} \right)^{2}},} & (5) \end{matrix}$

and τ, as above, is the decay time constant of the accelerating structure 22.

As known to persons of ordinary skill in the art, the coupling coefficient of the accelerating structure 220 (also referred to as an “accelerating resonator,” or a “resonator,” or “RF accelerating structure”) to the external circuit (feeding waveguide 222) depends on the current accelerated in (and interacting with) the resonator 220. Typically, the presence of current decreases the coupling coefficient that is measured through VSWR (voltage standing wave ratio), and the phase of the reflected signal from the resonator 220. Initially (without current), the resonator 22 needs to be over-coupled and to have a coupling coefficient of greater than β=1. The optimum coupling coefficient β₀ is a value that allows the resonator 22 to be matched with an external waveguide 222 at the accelerating current I. When the coupling coefficient β is equal to β₀, the coupling is referred to herein as “optimal.” An exact calculation of the optimum β₀ may be found by reference to Sobenin et al., Electrodynamic Characteristics of Accelerating Cavities (Eng. trans.), CRC Press, particularly at p. 121 (Eqn. 4.49), (1999), Collin, Foundations for Microwave Engineering, McGraw-Hill, (1^(st) ed., 1992), and Gao, Analytical formula for the coupling coefficient β of a cavity waveguide coupling system, Physics Research A, vol. 309, pp. 5-10 (1991), all of which are incorporated herein by reference.

SUMMARY OF EMBODIMENTS OF THE INVENTION

In accordance with embodiments of the present invention, a method is provided for generating a multiple-energy X-ray pulse. The method has steps of:

-   -   a. generating a beam of electrons with an electron gun;     -   b. modulating the beam of electrons prior to injection into an         accelerating structure to achieve at least a first specified         beam current amplitude and a first specified beam current         temporal profile, and a second specified beam current amplitude         and a second specified beam current temporal profile, the beam         of electrons characterized by an electron beam pulse duration;     -   c. applying to the accelerating structure a radio frequency         field with a specified RF field amplitude and a specified RF         temporal profile characterized by an RF pulse duration;     -   d. injecting the beam of electrons with the first specified beam         current amplitude and then with the second specified beam         current amplitude after a specified delay, in such a manner as         to achieve at least two distinct endpoint energies of electrons         accelerated within the accelerating structure during a course of         a single RF-pulse;     -   e. optimizing a coupling coefficient between an RF source and         the accelerating structure so as to achieve zero RF power         reflection at the specified beam current;     -   f. accelerating the beam of electrons with the radio frequency         field within the accelerating structure to produce accelerated         electrons; and     -   g. impinging the accelerated electrons upon a target for         generating Bremsstrahlung X-rays.

In accordance with alternate embodiments of the present invention, the electron beam pulse duration may be shorter than the RF pulse duration by a sum of onset delays defined by a filling time of the accelerating structure. The first specified beam current temporal profile and the specified RF temporal profile may begin substantially contemporaneously. The beam of electrons may be characterized by an amplitude of injection current, the amplitude of which may decrease during the multiple-energy X-ray pulse.

In accordance with further embodiments of the present invention, an end-point energy characterizing the beam of electrons may increase during the multiple-energy X-ray pulse. Applying the RF field may include modulating an RF source, such as by varying at least one of an input voltage and an input current to the RF source. Applying the RF field may also include modulating an RF modulator disposed between an RF source and the accelerating structure while the RF-source provides a constant level of power at constant frequency. A high-energy portion of the multiple-energy X-ray pulse may be characterized by a lower electron beam flux than a low-energy portion of the multiple-energy X-ray pulse.

In accordance with yet further embodiments of the present invention, the accelerating structure may include a standing wave resonator. The RF source may provide a constant level of RF-power; and a temporal profile of the electron gun may be characterized by at least two distinct levels of the amplitude of electron beam I_(n), wherein each amplitude of the at least two distinct levels of amplitude is created in temporally descending order and wherein the coupling coefficient of accelerating resonator β₀ is chosen to be optimal at a first level of beam current I₁. Breaking points, as defined below, pertaining to at least two distinct levels may include dynamically variable breaking points.

In accordance with other embodiments of the present invention, a first portion of an electron beam current pulse may start with an optimum delay (t_(b1)) relative to a to beginning of an RF pulse and end at a specified time (t₁), and each successive portion of the electron beam current pulse may be characterized by a distinct level I_(n) and start with a

delay

$t_{bn} = {\tau \cdot {\ln \left( \frac{I_{n - 1}}{I_{n}} \right)}}$

corresponding to an end of a previous pulse with τ representing a decay time of the accelerator structure.

In yet other embodiments of the invention, an end point of each pulse t_(n) may include a dynamically variable end point. There may be additional steps of

-   -   a. supplying an RF-power pulse to the accelerating structure at         two distinct energy levels, P_(L), P_(H), corresponding,         respectively, to low and high energy portions of the beam of         electrons;     -   wherein an electron gun temporal profile is characterized by two         distinct levels of the amplitude of electron beam, I_(L), I_(H),         corresponding, respectively, to beam current levels for         achieving low and high energy portion of the beam of electrons;         and     -   b. optimizing a coupling coefficient of accelerating resonator         β₀ based upon parameters of the low energy beam.

One additional embodiment of the invention has further steps of:

-   -   a. forming an ascending order of two distinct X-ray energy         levels by applying a corresponding ascending sequence of         RF-power levels and a descending sequence of injection current         pulses;     -   b. delaying a low energy component of the pulse by a low energy         current delay t_(bL) defined by

${t_{bL} = {\tau \cdot {\ln \left( \frac{2 \cdot \beta_{0}}{\beta_{0} - 1} \right)}}};$

and

-   -   c. delaying a high energy component of the pulse by a high         energy current delay t_(bH α)defined by

$t_{{bH}\; \_ \; a} = {\tau \cdot {\left\lbrack {{\ln \left( \frac{I_{L}}{I_{H}} \right)} + {\ln\left( {1 + {\frac{2 \cdot \beta_{0}}{\beta_{0} - 1}\sqrt{\frac{P_{H}}{P_{L}} - 1}}} \right)}} \right\rbrack.}}$

Another additional embodiment of the invention has further steps of:

-   -   a. forming a descending order of two energy levels by applying a         corresponding descending sequence of RF-power levels and an         ascending sequence of injection current pulses;     -   b. delaying a high energy component of the pulse by a high         energy current delay t_(bH) defined by equation

${t_{bH} = {\tau \cdot {\ln\left( \frac{\sqrt{4\beta_{0}{rLP}_{H}}}{I_{H}{rL}} \right)}}};$

and

-   -   c. delaying a low energy component of the pulse current delay         t_(bL_d) is defined by

$t_{{bL}\; \_ \; d} = {\tau \cdot {\left\lbrack {{\ln \left( \frac{I_{H}}{I_{L}} \right)} + {\ln\left( {1 + {\frac{\sqrt{4\beta_{0}{rLP}_{H}}}{I_{H}{rL}}\sqrt{1 - \frac{P_{L}}{P_{H}}}}} \right)}} \right\rbrack.}}$

BRIEF DESCRIPTION OF THE DRAWINGS

The foregoing features of the invention will be more readily understood by reference to the following detailed description, taken with reference to the accompanying drawings, in which:

FIG. 1 depicts a typical high-energy transmission X-ray inspection system, in the context of which embodiments of the present invention are advantageously applied.

FIG. 2 shows a block diagram of an X-ray source employing an accelerating structure and modulated current injection and RF excitation. A dual-energy pulse is created by applying the I_(L) and I_(H) currents to two distinct portion of the single RF pulse.

FIG. 3 depicts dependence of energy and current within a microwave pulse, in accordance with the same or another embodiment of the present invention.

FIG. 4 shows a block diagram of an X-ray source employing a standing wave accelerating structure with modulated current injection and RF excitation, in accordance with the present invention.

FIG. 5 shows a linac implementation with a dual-energy pulse created by applying the I_(L) and I_(H) currents with optimal delays, in accordance with an embodiment of the present invention.

FIG. 6 shows a block-diagram of an intra-pulse dual-energy linac where the energy modulation is achieved by varying both input RF-power and injection current provided into a standing wave accelerating system, in accordance with an embodiment of the present invention.

FIG. 7 shows an example of creating an ascending order of beam energy steps W_(n) by applying an ascending order of RF-power levels P_(n) with corresponding descending order of injection currents I_(n), in accordance with an embodiment of the present invention.

DETAILED DESCRIPTION OF EMBODIMENTS OF THE INVENTION

Definitions. The term “multiple-energy” shall refer to an X-ray inspection modality in which differential transmission through a medium by X-rays of distinct spectral composition is used to characterize the medium.

The term “pulse duration,” denoted t_(RF), refers to the duration of time that RF excitation is applied to a linac accelerating structure.

The term “breaking point” of a current pulse of duration t_(L)+t_(H) has a very specialized meaning herein: It is defined to be the value of t_(L)/t_(H), where t_(L) and t_(H) refer, respectively, to durations of the current pulse during which an electron beam emitted from an accelerating structure is characterized by low- and high-energies, respectively.

In the case of a current pulse having multiple sub-pulses, any pair of sub-pulses may be characterized by a breaking point that is defined as the value of t_(L)/t_(H), where t_(L) and t_(H) refer, respectively, to durations of current sub-pulses.

A breaking point of a set of current pulses or sub-pulses shall be said to be “dynamically variable” if the breaking point is adapted to be varied from one pair of current pulses or sub-pulses to another pair of current pulses or sub-pulses. Similarly, and end point of a pulse may also be characterized as “dynamically variable” if the end point is adapted to be varied from one current pulse to another current pulse.

The term “current level,” as it pertains to an x-ray source, refers to an average flux of electrons incident upon the target, expressed in milliAmperes (mA), and averaged over a specified duration of time. Unless otherwise indicated, the specified duration over which an average is taken is the duration of a pulse.

The term “current amplitude,” as it pertains to an x-ray source, refers to a value of an instantaneous flux of electrons incident upon the target, expressed in mA.

The term “onset delay” refers to a period between application of an RF field to an accelerating structure and injection of a pulse or sub-pulse of electron current into the accelerating structure. Where there are multiple sub-pulses of electron current, then the cumulative onset delays of the respective sub-pulses make up a sum of onset delays.

As used herein, the term “Bremsstrahlung” may be used to denote X-ray emission produced through impingement of high-energy electrons onto a metallic target, and, additionally, the physical process involved in that phenomenology.

The systems and methods described herein may be described in terms of X-rays, however the applicability of the teachings to other spectral ranges is clear, and encompasses, within the scope of the invention, all manner of penetrating radiation.

Various embodiments of the invention described herein employ variation of the spectral content of an X-ray pulse during the course of the pulse to discriminate differences in X-ray transmission of a medium in different energy regimes. Approaches taught in accordance with the present invention are particularly advantageous in cases where high speed of scanning is required, such as train or high-throughput scanners.

In accordance with embodiments of the present invention, a novel apparatus and novel methods are provided that may advantageously create a pulse profile of a multi-energy beam in such a manner that would improve material discrimination while preserving the highest possible scanning speed and allow optimizing dose to cargo and environment.

A novel source of penetrating radiation, designated generally by numeral 20, in accordance with an embodiment of the present invention, is now described with reference to FIG. 2. A linac 21, depicted in FIG. 2, includes accelerating structure 22, and an electron gun 23 serving as an injector of electrons emitted by cathode 235. Used in conjunction with linac 21 are an X-ray target 24, an RF-source 25, an RF-source modulator 26, an RF isolator 27 and an electron gun modulator 28. An RF-circuit 29 provides a constant level of microwave power within pulse duration t_(RF). Electron gun 23, driven by e-gun modulator 28, provides an electron beam 220 characterized by a two-level injection current pulse 210 (also referred to herein as “the pulse”) into accelerating system 20 with a total duration t_(p)≤t_(RF). Injection current pulse 210 may also be referred to herein as “injection current,” and its amplitude, which, as defined above, corresponds to its instantaneous value of electron flux in mA, denoted I_(L).

The first portion 212 of the pulse 210, characterized by a higher amplitude of the injection current I_(L) due to higher beam loading, creates a low energy portion of the beam pulse, where the low energy portion is designated by W_(L). (For avoidance of ambiguity, it is to be noted that W_(L) refers both to the low energy portion of the pulse, and to the value of the instantaneous endpoint energy characterizing the low energy portion of the pulse. The same applies, mutatis mutandis, to W_(H), the high energy portion of the pulse.) The second portion 214 of the injection current pulse with lower amplitude I_(H) produces a high energy portion of the beam pulse W_(H). The “breaking point” of the pulse, which, as defined above, has the specialized meaning of the value of t_(L)/t_(H), may be variable, within the scope of the present invention, thereby enabling dynamic control of the dose of emitted X-rays to cargo and environment. In accordance with certain embodiments of the present invention, the breaking point may advantageously be varied from pulse to pulse, rendering it a dynamically variable breaking point.

Reference is made now to FIG. 3, where dependence of energy and current within a microwave pulse is depicted. Dashed line 32 shows the injection current, while solid line 34 represents the beam energy. Dotted line 30 shows energy dependence at constant current I_(H), as previously discussed. As the result of different beam loading effect in the front and rear portion of the pulse, the beam has two distinct energy levels, and, in one embodiment: W_(L(0.-1.5)μs)=3.9 MeV and W_(H (1.8-3.3) μs)=5.8 MeV. As used herein, “energy level” refers to the instantaneous end-point energy of an ensemble of photons, distributed in energy according to an essentially Bremsstrahlung spectrum of photon energies, or otherwise.

Optimizing coupling coefficient β₀ of the accelerating structure for the parameters of the single energy beam is known in the prior art, and has been described in the Background Section above.

In FIG. 4, a linac 21 is shown that is similar to the one shown in FIG. 2. The accelerating system 20 is based on a standing wave structure 42 (otherwise referred to herein as a “standing wave resonator”), the distinction of which with respect to a traveling wave structure has been laid out by Miller (1986). The coupling coefficient β₀ is chosen to be optimal at the current I_(L), using the algorithm that was laid out in detail above. The value of I_(L) is chosen to provide energy W_(L), and this energy value W_(L) remains constant over an entire sub-pulse duration if I_(L) is applied with delay t_(bL) with respect to the beginning of the RF pulse. Low energy current delay t_(bL)is defined by Eq. (4) above. At the end of the low energy pulse (after t_(L)), the current turns off. The value of I_(H) is chosen to provide W_(H), and this energy level remains constant if I_(H) current is applied with delay t_(bH) counting from the end of low energy pulse t_(L). The high energy pulse delay is defined by equation:

$\begin{matrix} {t_{bH} = {\tau \cdot {{\ln \left( \frac{I_{L}}{I_{H}} \right)}.}}} & (6) \end{matrix}$

In accordance with certain embodiments of the present invention, the t_(L) point may be allowed to vary, thereby allowing the ratio t_(L)/t_(H) (defined herein as the “breaking point”) to be varied, and thus advantageously providing for dynamic control of the X-ray dose to cargo and environment.

The average current during the lower energy portion of the pulse will be referred to herein as the low energy current, and, mutatis mutandis, the average current during the higher energy portion of the pulse will be referred to herein as the high energy current.

The energies within each portion of the pulse will remain constant as long as the “low energy” current begins to be applied after a delay of t_(bL), and as long as the “high energy” current is applied with a delay of t_(bH). The constancy of energy within each of the LE and HE portions of the pulse is beneficial for material discrimination: the energy spectrum of X-ray beam remains constant hence no additional calibration point(s) is required.

A multi-energy pulse configuration, with greater than two distinct energies during the duration of each pulse, may be created in a similar fashion to that described above, using a standing wave accelerating structure.

-   -   The coupling coefficient β₀ is chosen to be optimal at current         I₁, using the design algorithm described in detail above;     -   The current I₁ has the highest value in the sequence, it is         applied first with delay of t_(b1) (defined by Eqn. (4)) thus         creating the lowest energy of the beam sequence.     -   The values of next current pulses I_(n) are preferably created         in descending order, thereby providing an ascending order of the         beam energy levels.     -   In preferred embodiments of the invention, after each sub-pulse         the current is turned off.     -   The current I_(n) for n-th sub-pulse is applied with delay         t_(bn), given by:

$\begin{matrix} {t_{bn} = {\tau \cdot {{\ln \left( \frac{I_{n - 1}}{I_{n}} \right)}.}}} & (7) \end{matrix}$

Each energy level end point t_(n) may still be allowed to vary, thereby advantageously providing for dynamic control of the dose to cargo and environment. The energy within each portion will remain constant as long as the lowest I₁ current is applied after a delay of t_(b1) (defined by Eqn. (4)) and each successive current step is applied with a delay of t_(bn).

An example of linac implementation with a dual-energy pulse that is created by applying the I_(L) and I_(H) currents with optimal delays is shown in FIG. 5. Parameters of the linac are identical to those that have been shown as an example in FIG. 3. Low energy current 51 is applied at an optimal delay of t_(bL)≅0.34 μs and is turned off after t_(L)≅1.45 μs. High energy current 52 is applied with delay of t_(bH)≅0.36 μs at 1.8 μs and is turned off at the end of microwave pulse.

FIG. 6 depicts a block-diagram of an intra-pulse dual-energy linac 60 where energy modulation is achieved by varying both input RF-power P_(H) and P_(L), as well as injection current I_(H) and I_(L) provided into the standing wave accelerating structure 42. Low energy current 61 and high energy current 62 are applied with optimum delays. RF-power modulation can be achieved by several known methods such as varying input voltage and current of RF-source 25, or varying the input RF-power to the RF-source by means of RF-source modulator 26, or else by manipulating output power of the RF-source with switches or regulators, generically referred to herein as RF modulator 255.

In accordance with another embodiment of the present invention, an ascending order of energy levels is created with constant amplitudes. As used herein, the term “amplitude” refers to an instantaneous flux of electrons within electron beam. In this embodiment.

-   -   A corresponding ascending sequence of RF-power levels and         descending sequence of injection current pulses is applied.     -   The coupling coefficient β₀ is chosen to be optimal at current         I_(L), using the design algorithm described in detail above.     -   A low energy current delay t_(bL) is as defined by Eqn. (4).     -   A high energy current delay t_(bH_a) is as determined using

$\begin{matrix} {t_{{bH}\; \_ \; a} = {\tau \cdot {\left\lbrack {{\ln \left( \frac{I_{L}}{I_{H}} \right)} + {\ln\left( {1 + {\frac{2 \cdot \beta_{0}}{\beta_{0} - 1}\sqrt{\frac{P_{H}}{P_{L}} - 1}}} \right)}} \right\rbrack.}}} & (8) \end{matrix}$

P_(H) and P_(L) refer, respectively, to RF power applied during high energy and low energy portions of the beam current pulse.

While using power modulation, certain benefits may be achieved by creating a descending order of energy levels with constant amplitudes. For implementing such the descending sequence:

-   -   A corresponding descending sequence of RF-power levels and         ascending sequence of injection current pulses is applied.     -   β₀ is chosen to match the accelerating system with low energy         current using design algorithm described in detail above.     -   A high energy current delay t_(bH) is determined using Eqn. (3)         using I_(H) and P_(H).     -   A low energy current delay t_(bL_d) is determined using

$\begin{matrix} {t_{{{bL}\; \_ \; d}\;} = {\tau \cdot \left\lbrack {{\ln \left( \frac{I_{H}}{I_{L}} \right)} + {\ln\left( {1 + {\frac{\sqrt{4\beta_{0}{rLP}_{H}}}{I_{H}{rL}}\sqrt{1 - \frac{P_{L}}{P_{H}}}}} \right)}} \right\rbrack}} & (9) \end{matrix}$

The block-diagram of the linac 71 depicted in FIG. 7 demonstrates the option of a linac with multi-energy pulses that is created by modulating both input RF-power and injection current. Linac 71 includes similar subsystems to those presented in FIGS. 2 and 4. FIG. 7 shows an example of creating an ascending order of beam energy steps W_(n) by applying an ascending order of RF-power levels P_(n)with a corresponding descending order of injection currents I_(n).

Embodiments of the invention described above are intended to be merely exemplary; numerous variations and modifications will be apparent to those skilled in the art. All such variations and modifications are intended to be within the scope of the present invention as defined in any appended claims. 

I claim:
 1. A method for generating a multiple-energy X-ray pulse, the method comprising: a. generating a beam of electrons with an electron gun; b. modulating the beam of electrons prior to injection into an accelerating structure to achieve at least a first specified beam current amplitude and a first specified beam current temporal profile, and a second specified beam current amplitude and a second specified beam current temporal profile, the beam of electrons characterized by an electron beam pulse duration; c. applying to the accelerating structure a radio frequency field with a specified RF field amplitude and a specified RF temporal profile characterized by an RF pulse duration; d. injecting the beam of electrons at a first specified beam current amplitude and then at the second specified beam current amplitude after a specified delay, in such a manner as to achieve at least two distinct endpoint energies of electrons accelerated within the accelerating structure during a course of a single RF-pulse; e. optimizing a coupling coefficient between an RF source and the accelerating structure so as to achieve zero RF power reflection at the specified beam current; f. accelerating the beam of electrons with the radio frequency field within the accelerating structure to produce accelerated electrons; and g. impinging the accelerated electrons upon a target for generating X-rays by Bremsstrahlung.
 2. A method in accordance with claim 1, wherein the electron beam pulse duration is shorter than the RF pulse duration by a sum of onset delays defined by a filling time of the accelerating structure.
 3. A method in accordance with claim 1, wherein the first specified beam current temporal profile and the specified RF temporal profile begin substantially contemporaneously.
 4. A method in accordance with claim 1, wherein the a beam of electrons is characterized by an amplitude of injection current, and wherein the amplitude of injection current decreases during the multiple-energy X-ray pulse.
 5. A method in accordance with claim 1, wherein an end-point energy characterizing the beam of electrons increases during the multiple-energy X-ray pulse.
 6. A method in accordance with claim 1, wherein applying the RF field includes modulating an RF source.
 7. A method in accordance with claim 6, wherein modulating the RF source includes varying at least one of an input voltage and an input current to the RF source.
 8. A method in accordance with claim 1, wherein applying the RF field includes modulating an RF modulator disposed between an RF source and the accelerating structure while the RF-source provides a constant level of power at a constant frequency.
 9. A method in accordance with claim 1, wherein a high-energy portion of the multiple-energy x-ray pulse is characterized by a lower electron beam flux than a low-energy portion of the multiple-energy X-ray pulse.
 10. A method in accordance with claim 1, wherein the accelerating structure includes a standing wave resonator.
 11. A method in accordance with claim 10, wherein the RF source provides a constant level of RF-power; wherein a temporal profile of the electron gun is characterized by at least two distinct levels of the amplitude of electron beam I_(n); wherein each amplitude of the at least two distinct levels of amplitude is created in temporally descending order; and wherein the coupling coefficient of accelerating resonator β₀ is chosen to be optimal at a first level of beam current I₁.
 12. A method in accordance with claim 11, wherein breaking points of said at least two distinct levels include breaking points that are dynamically variable.
 13. A method in accordance with claim 11, wherein a first portion of an electron beam current pulse starts with an optimum delay (t_(b1)) relative to a beginning of an RF pulse; wherein said electron beam current pulse ends at a specified time (t₁); and wherein each successive portion of the electron beam current pulse is characterized by a distinct level I_(n) and starts with a delay $t_{bn} = {\tau \cdot {\ln \left( \frac{I_{n - 1}}{I_{n}} \right)}}$ corresponding to an end of a previous pulse t_(n-1), with representing a decay time of the accelerator structure.
 14. A method in accordance with claim 12, wherein an end point of each pulse t_(n) includes a dynamically variable end point.
 15. A method in accordance with claim 11, further comprising: a. supplying an RF-power pulse to the accelerating structure at two distinct energy levels, P_(L), P_(H), corresponding, respectively, to low and high energy portions of the beam of electrons; wherein an electron gun temporal profile is characterized by two distinct levels of the amplitude of electron beam, I_(L), I_(H), corresponding, respectively, to beam current levels for achieving low and high energy portion of the beam of electrons; and b. optimizing a coupling coefficient of accelerating resonator β₀ based upon parameters of the low energy beam.
 16. A method in accordance with claim 15, further comprising: c. forming an ascending order of two distinct X-ray energy levels by applying a corresponding ascending sequence of RF-power levels and a descending sequence of injection current pulses; d. delaying a low energy component of the pulse by a low energy current delay t_(bL) defined by ${t_{bL} = {\tau \cdot {\ln \left( \frac{2 \cdot \beta_{0}}{\beta_{0} - 1} \right)}}};$ and e. delaying a high energy component of the pulse by a high energy current delay t_(bH a) defined by $t_{{bH}\; \_ \; a} = {\tau \cdot {\left\lbrack {{\ln \left( \frac{I_{L}}{I_{H}} \right)} + {\ln\left( {1 + {\frac{2 \cdot \beta_{0}}{\beta_{0} - 1}\sqrt{\frac{P_{H}}{P_{L}} - 1}}} \right)}} \right\rbrack.}}$
 17. A method in accordance with claim 15, further comprising: c. forming a descending order of two energy levels by applying a corresponding descending sequence of RF-power levels and an ascending sequence of injection current pulses; d. delaying a high energy component of the pulse by a high energy current delay t_(bH) defined by equation ${t_{bH} = {\tau \cdot {\ln\left( \frac{\sqrt{4\beta_{0}{rLP}_{H}}}{I_{H}{rL}} \right)}}};$ and e. delaying a low energy component of the pulse current delay t_(bL_d) is defined by $t_{{bL}\; \_ \; d} = {\tau \cdot {\left\lbrack {{\ln \left( \frac{I_{H}}{I_{L}} \right)} + {\ln\left( {1 + {\frac{\sqrt{4\beta_{0}{rLP}_{H}}}{I_{H}{rL}}\sqrt{1 - \frac{P_{L}}{P_{H}}}}} \right)}} \right\rbrack.}}$ 